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A numerical analysis shows that a class of scalar-tensor theories of gravity with a scalar field 
minimally and nonminimally coupled to the curvature allows static and spherically symmetric black 
hole solutions with scalar-field hair in asymptotically flat spacetimes. In the limit when the horizon 
radius of the black hole tends to zero, regular scalar solitons are found. The asymptotically flat solu- 
. tions are obtained provided that the scalar potential V(<f>) of the theory is not positive semidefinite 

' and such that its local minimum is also a zero of the potential, the scalar field settling asymptoti- 

| cally at that minimum. The configurations for the minimal coupling case, although unstable under 

spherically symmetric linear perturbations, are regular and thus can serve as counterexamples to 
the no-scalar-hair conjecture. For the nonminimal coupling case, the stability will be analyzed in a 
forthcoming paper. 

. PACS numbers: 04.70.Bw, 04.20.Jb, 97.60.Lf 

It was more than thirty years ago when Ruffini and Wheeler l) proposed the so called no-hair conjecture for black 
■ holes. This conjecture states that black holes are completely characterized by their mass, charge, and angular mo- 
mentum. In order to prove this conjecture, some people established several no-hair theorems in theories which couple 
classical fields to Einstein gravity, notably the black hole uniqueness theorems in Einstein-Maxwell (EM) theory 
which establishes that all the solutions of black holes in EM theory are stationary and axially symmetric and contained 
within the Kerr-Newman family. Other theorems proved by Chase Bekenstein Hartle [j| and Teitelboim [||, 
showed that stationary black holes solutions are hairless in a variety of theories coupling different classical fields to 
Einstein gravity. A key ingredient in the no-hair-theorem proofs relies on the assumption of asymptotic flatness (AF) 
and on the nature of the energy-momentum tensor. For instance, in recent years counterexamples to the no-hair 
conjecture were found in several theories with non-Abelian gauge fields which include the Einstein- Yang-Mills Q, 
'. Einstein- Yang-Mills-Higgs (EYMH) @], Einstein- Yang-Mills dilaton (EYMD) @, the Einstein-Skyrme [lj| and the 
jpP" Einstein-non-Abelian-Procca Q theories. The existence of solitons and hairy black holes in these theories i s 
J> , associated with the non-Abelian gauge fields present in them I n the case of matter composed by a single scalar 
field minimally coupled (MC) to gravity, Sudarsky |l3j |. has proved a very simple no- hair theorem. Among other as- 
sumptions, the validity of that theorem relies on the fulfillment of the weak energy condition (WEC) which constrains 
the scalar potential to be positive semi-definite. Recently a more general no-hair theorem for black holes has been 
proved. This rules out a multicomponent scalar field coupled minimally to gravity satisfying the WEC, but its field 
Lagrangian is not quadratic in the field derivatives |14| . In fact a central argument for probing all these theorems 
relies on the fulfillment of the WEC which is generally believed to be satisfied for all physically reasonable classical 
matter 

In the case of nonminimally coupled (NMC) scalar fields, there are only a few no-hair theorems [T(l Il7j. The AF 
hairy black hole solution corresponding to a scalar field conformally coupled to gravity (the Bronnikov-Melnikov- 
Bocharova-Bcckenstein solution [l8|) which is often used as a counterexample for the failure of the no-scalar-hair 
conjecture, has been criticized as a non genuine black hole solution in that the Einstein field equations are not verified 
at the horizon |l9j, and moreover, in that it is not a solution at all if one demands a bounded scalar field throughout 
the static region [12(a)]. 

Recently, in the context of asymptotically AdS spacetimes, scalar-hairy black holes (SHBH) were found for some 
scalar- field potentials V{4>) J2jj. A further analysis along this line [2lJ, suggested in addition that by adjusting the 
parameters of V(<j>), the effective cosmological constant could cancel out so that asymptotically flat SHBH could exist. 
Encouraged by the findings in the AdS context, we show in this letter numerical evidence of SHBH solutions in 
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static, spherically symmetric and asymptotically flat spacetimes. The solutions are regular throughout the static 
region (from the Killing horizon to spatial infinity). In particular, for the MC case, the main assumption for the 
Sudarsky's no-hair theorem 0] not to be valid here is that the class of scalar-field potentials V((j>) used to construct 
the solutions are not positive-semidefinite (the WEC is violated). Moreover, the class V(<f>) is in fact such that 
V(4>c) = = d c j,V\ t p c and d' 2 J( pV\^ c > 0, that is, a root of V{<f>) and a local minimum are located at the same place. 
This feature and the fact that the non-trivial scalar field settles asymptotically at the local minimum (and therefore 
also at a root) of V ((/)), leads to configurations that truly represent asymptotically flat solutions. When taking the 
limit r/j — * (r/j being the horizon radius), the SHBH tend to regular scalar solitons (scalarons). Such results can 
extend to the NMC case as it turns numerically (the no-hair theorems for the NMC are avoided by considering the 
scalar potentials as mentioned above). 

The SHBH solutions although genuine asymptotically flat solutions, and perhaps the first examples with a regular 
horizon and with an explicit scalar field potential, constitute only a "weak" counterexample to the no-scalar-field-hair 
conjecture since they turn to be unstable with respect to radial perturbations. 

Einstein- scalar field equations and Lagrangian.- We will consider a model of a scalar field NMC to gravity and with 
a potential. The simplest model of this kind is obtained by considering the Lagrangian 



L = 



-9 



lD7r 2 



V{4>) 



(i) 



By choosing the NMC constant £ = 0, the above theory corresponds to the minimal coupling case, and £ = —1/12 
corresponds to the conformally coupled scalar field (units where Go = c = 1 are employed). The gravitational field 
and scalar field equations following from the Lagrangian Q can be written as 



= 8nT^ , Of, + 2£R = 



dV(4>) 



(2) 



where 



T MV = G e ff \ (V^V^ - M „ 



4£ [V M (0V v 0) - g^VxiW* 



(3) 



is an effective energy-momentum tensor which includes all the contributions of the scalar field, and G c s is an effective 
gravitational "constant" which explicitly depends on the scalar field: 



G 



1 



cff 



(1 + 16tt£0 2 ) ■ 

We will focus on a metric describing spherical and static spacetimes: 

ds 2 = -Ne 25 dt 2 + N~ x dr 2 + r 2 (d8 2 + sin 2 



(4) 



(5) 



where N = (1 — 2m(r)/r) and the functions m and S depend only on the coordinate r. For the scalar field we also 
assume <p — 4>{r). The resulting field equations are 



d r m = Anr 2 E , 

4:7TV 

d r 6 = —[E + S r 



d 



dV{<t>) 



d4> 



2&R 



(6) 
(7) 

(8) 



where E = — T/, S r = T T r and R are the local energy density, the radial pressure and the Ricci scalar, respectively, 
given by the following expressions containing no-second order derivatives: 

G c ff 



E = 

S r = 
R = 



1 + 167r£r</>(d r 0)G eff 

G e ff 

1 + 1927r£ 2 (/> 2 Gcff 
G c ff 



= 167r£r0(d r 0)G c ff V[<p) 



N(d r t) 2 



1 + 8£ + 640 J G e ff) + 4£0 



V 2 

dV(4>) 



V(<t>) (1 - 64< 2 ^ 2 G eff ) 



1 + 16Tr^rcj)(d r (j))G eS 
8ttG c s 



-Atj<t>(d r <j>)m N{d r cj>) 2 S^jdr^N 

12 + — o v w : 



1 + 192< 2 2 G eff 



N{d r 4>y{l + 120 + 4V(4>) + 12£0 



dV{/f>) 



(9) 
(10) 

(11) 
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Boundary conditions and numerical methodology.- For black hole configurations we demand regularity on the event 
horizon r*/, and for scalarons regularity at the origin r = 0. This implies the following conditions for the fields 
(hereafter, the subscripts 7i' stand for a quantity to be evaluated at r^): 

m h = — , 6{r h ) = S h , <j){r h ) = <j> h , {d r <j)) h = — — f— , (12) 



where 



^rl(S r h -E h ) + l 



Rh - i^^m [Vh+3 ^ v)h] - (13) 



32ttG^ 



S > 1 ~ Eh = ~ l + 192rtjW? h [^{d^V) h + V h (1 - 6^e^G h cS )} - G h cS V h . (14) 



The value Sh is fixed so that the desired asymptotic behavior (see below) is obtained. For the scalaron case, regularity 
at the origin r = results by taking the limit rh — > in the above regularity conditions. In addition to the regularity 
conditions, we impose asymptotically flat conditions on the space-time (for black holes and scalarons). These imply 
the following conditions on fields when r — > oo, 

ra(oo) = Madm , 8(oo) = , 4>(oo) = ^ . (15) 

Here Madm, is the ADM- mass associated with a given configuration. Actually, the value (poo will correspond to the 
local minimum (which is also a root) of V(4>) (see below). 

For a given scalar field potential and for a fixed £ (i.e., for a given theory), the family of SHBH configurations will be 
parametrized by the arbitrary free parameter r/j which specifies the location of the black hole horizon. As explained 
below, (f>h is a shooting parameter rather than an arbitrary boundary value, and its value is determined so that the 
above AF conditions are fulfilled. Therefore for SHBH, Madm = Madm^/i)- As mentioned above, the scalarons are 
contained as a limiting case when r/,-*0 (in this limit, 4>q, the value of </> ar r — 0, is the shooting parameter), and 
the corresponding configuration is characterized by a unique Madm- 

Using the above system of field equations together with the regularity and asymptotic conditions, we have performed 
a numerical analysis for one class of scalar field potential and for different values of £. 

Numerical results. We choose the following asymmetric scalar-field potential leading to the desired asymptotically 
flat solutions: 



{(p-a) (0 - a) + 2771772 



{cp-af , (16) 



where A, r\i and a are constants. For this class of potentials one can see that, for 771 > 2772 > 0, 4> = a corresponds to a 
local minimum, cf> — a + r/i is the global minimum and <fi — + 772 is a local maximum. The key point in the shape of the 
potential, V(4>), for the asymptotically flat solutions to exist, is that the local minimum = V(a) is also a zero of 
V(4>) |22^ . The dynamics of scalar fields with such a potential has been analyzed in the past within the aim of stud ying 
quantum tunneling from the false vacuum through the true vacuum (see Ref. 23] ) . A particle- mechanics analogy [23| , 
helps to understand the existence of classical solutions that interpolate between the local minimum <f> = a and a value, 
4> = 0/1.0 , near the global minimum, and rolling across the local maximum. This requires a suitable shooting method 
in order to find the correct value cj>h f° r the scalar field to reach the local minimum asymptotically j24[. We have 
enforced such a method here. The asymptotic behavior of the scalar field shows that <p(r) ~ a + const. e ±%/Al ' 1 ' )2 ''/r. 
The shooting method allows to eliminate the runaway solutions. 

Since V((f>) is not positive definite (we assume A > 0), the WEC is violated. On the other hand, taking potentials 
with extrema different from zero (V^-£ 7^ 0), leads to solutions that are asymptotically AdS (if V™^ < 0), or 
asymptotically de Sitter (if > and provided that suitable boundary conditions are imposed at the cosmological 
event horizon) with acting as an effective cosmological constant A c fj (cf. Ref. [2(J for an analysis with A e g 7^ 0). 

For the numerical analysis, A determines the scale of different quantities. In particular r c = l/\/A has been used 
as a lenght scale. 

Case £ = 0. Figures ^ an d 121 show examples of numerical asymptotically flat SHBH and scalaron solutions 
respectively, for the potential V(<j>). As depicted by those figures, a generic behavior of the solution for <f> is that 
it decays exponentially (0 ~ const. e _v ' A '' 1 ' ,2r /r) to the local minimum Vj]^ before reaching the asymptotic value at 
4> = a (a = in these examples). <ph,o is the shooting value for which <f> rolls up to <f> — a asymptotically. The right 
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FIG. 1: 

Black hole configuration constructed with £ = and V((f>) as given by Eq. 1161 with parameters 771 = 0.5, 772 = 0.1, a — 0, and 
rh — 0.1/\/A, <f>h ~ 0.40786. The left and middle panels depict the scalar field and the mass function respectively. The latter 
converges to A/adm ~ 3.843/vA. The right panel depicts the metric potentials \J—gtt (solid line), yjg r r (dashed line), e 5 

(dash-dotted line) and 5 (dotted line). 
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FIG. 2: 

Same as Fig^for the soliton case. Here <j>o ~ 0.40594 and r c = 1/vA. Note that the solutions are globally regular, notably at 

the origin, and Madm ~ 3.827/ V\. 



panel of ngQ] (figEJ) shows that the SHBH (scalaron) solution asymptotically approximates the Schwarzschild BH 
(Schwarzschild "exterior" ) solution with mass A/adm respectively. 

The ADM mass, Madm('"/i)) f° r different values of the parameters (A, 771,772), turns to be larger than the corre- 
sponding mass of the Schwarzschild BH while keeping fixed ru ■ It seems therefore that the lowest bound for Madm is 
M^dm which corresponds to the hairless Schwarzschild BH. In the limit — > 0, the lowest bound is then Madm = 
which is no other but the trivial soliton of the Minkowski spacetime. The theory also admits the no-hairy solutions 
4>(r) — </>^°j al (4>{r) — 0max)i f° r which tf> settles at the global minimum (local maximum) ^ 1 i ° bal (V^jf 1 ) respec- 
tively. However, those solutions are not asymptotically flat, but rather correspond to the Schwarzschild AdS solution 
(Schwarzschild dS solution) with V^f;° bal (Vj]^ 1 ) playing the role of a negative (positive) cosmological constant. 

Since Mp^m{ r h) > ^ADu( r h)^ i- e -j Madm^) is bounded from below (while keeping fixed Th which in these 
coordinates is equivalent to fixing the area, Ah, of the horizon) with M^m corresponding to the Schwarzschild black 
hole, heuristically one expects that the hairy BH configurations found here with fixed AF boundary conditions are 
unstable, since the stable ones would correspond to those with the lowest Madm^/i) (cf. Ref. H^]). Such a behavior is 
a signature of the unstable nature of the SHBH configurations which is confirmed by a perturbation analysis (similar 
arguments can be applied to the scalarons to exhibit their unstable nature). 

Case £ 7^ 0. After performing numerous experiments for the minimal coupling case £ = 0, we also found (for 
some set of values £ > 0) that asymptotically flat hairy black holes and scalarons exist for the NMC case as well (the 



5 



resulting figures are qualitatively similar to Figs ^ and [2J. However, a more exhaustive analysis is to be performed in 
this case. For instance, the stability in the NMC is an issue that remains to be investigated. This will be the object 
of a forthcoming paper in which we also plan to provide a full sample of black hole and soliton configurations, and to 
analyze more carefully the quantitative behavior of Madm = Mabm (fh > for a larger range of £ |27| . 

Linear stability analysis.- We now turn to the question of the stability of the asymptotically flat space-times 
presented here. For simplicity, the following analysis is limited to the case £ = 0. We consider only spherically 
symmetric time-dependent linear perturbations to the spherically symmetric static metric and the scalar field 
configuration discussed in the previous section. These perturbations will be parametrized by 

4>(t,r) = 0(r)+&(t,r) , (17) 
9tt(t,r) = g tt (r)(l-h (t,r)) , (18) 
g rr (t,r) = g rr (r)(l + h 2 (t,r)) , (19) 

where gtt{r) — —Ne 2S , g rr {r) = iV _1 , (f>(r) represent the non-perturbed solution and 4>i(t,r), h (t,r), /i 2 (i,r) denote 
a small perturbation to the non-perturbed solution. It is straightforward to obtain the metric perturbations as a 
function of the scalar field perturbation 

hi = 8irr(d r <j))(f>i , d r h = d r h 2 — l6Trr(d r (j))(d r (f>i) , (20) 

and the linear Schrodinger type equation for the scalar field perturbation 

dr p~ s 

-dl r ^ + v: s [nty = ,-£ = —, (21) 



where we have used ip = r<j>\, r*(r) is the "tortoise" coordinate, and the effective potential is given by 
V cS (r) = Ne 25 

Figure 13 depicts V e g for the SHBH (left panel) and soliton (right panel) respectively. 



!L Us + m - s M Us + r M. + l\ + Mt z + % 



(22) 




FIG. 3: 

The potential V4a corresponding to the SHBH solution of fig.l (left panel) and to the soliton solution of fig. 2 (right panel). 

One can seek mode perturbations ip(r) = x(r)e 4<Tt , so that the Equation (|2~T)> writes (— D a D a + V e s ) X = (j2 X 
where D a is the derivative operator associated with an auxiliary metric of a manifold M. According to the theorem 
proved in |26| , a sufficient condition for the static and spherically symmetric configurations to be unstable, is that the 
operator A = —D a D a + V e g be negative in the Hilbert space L 2 (M). This theorem can be easily implemented when 
the background (static) solution is given analytically. For instance, one can see explicitly, that the Schwarzschild 
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N 

2M 



1 — 2M^ DM /r is stable within these kind of perturbations and within this class of potentials, 
■^71 ^2 > in the region < r < 00. A result, which is heuristically confirmed by the 



0=0 

ADM 



solution S(r) = 0, 
since V c g(r) = N 

fact that A^adm ^ s a l° wer bound within the set of configurations with fixed |27| . However, when the background 
solution is given numerically the implementation of the theorem is less trivial since it requires the use of a suitable 
square integrable vector \f in L 2 (M) to compute < '5, >. Instead, we have opted to explicitly solve the associated 
Schrodinger-like equation of Eq. 12 111 with respect to the original r coordinate to find a "bound state" with negative a 2 . 
Imposing regularity at (origin), i.e., x|r h ,o = 0; an d demanding that x( r ) — * f° r r ~~ * 00 we found a 2 ~ —0.00241 
and a 2 ~ —0.00243 for the SHBH and scalarons solutions respectively (see Fig. 0J. This shows that the modes 
ip(r) — x(r)e ltTt can grow unboundedly with time, and therefore leading to the conclusion that the SHBH and soliton 
configurations are unstable. 



Asymptotically Flat Black Holes 
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FIG. 4: 

The mode x corresponding to the perturbation of the SHBH solution of fig.l (left panel) with a 2 ~ —0.00241 and the soliton 

solution of fig. 2 (right panel) with a 2 ~ —0.00243 respectively. 

Conclusions. We have found black hole solutions that support a nontrivial scalar hair in spherical, static and 
asymptotically flat spacetimes. Such solutions are perhaps the first regular ones throughout the static region (from 
the horizon to spatial infinity) with an explicit scalar field potential [28j. In the limit where — ► 0, regular scalar 
solitons (scalarons) are obtained. In the minimal coupling case £ = 0, the solutions turn to be unstable with respect to 
linear-spherical perturbations and their ADM mass turns to be larger than the corresponding Schwarzschild black hole 
while keeping fixed rj,. Therefore, these are to be regarded as weak counterexamples to the no-scalar-hair conjecture. 
For nonminimal couplings £ 7^ 0, similar solutions are found, however their stability remains to be investigated. 
Another analysis to be performed are the collapse of these hairy black holes violating WEC and the possibility of 
having hairy black holes with nodes in the scalar field using higher order potentials. An interesting outcome of the 
above results is that it opens the possibility for a further study of the solutions in the context of the isolated horizon 
formalism (cf. Ref. [29|). 

We are indebted to A. Corichi, D. Sudarsky and T. Zannias for discussions, and to E. Ayon for providing some 
references. We particularly thank O. Sarbach for his comments and suggestions on the stability analysis. U.N. would 
like to acknowledge partial support from SNI. M.S. acknowledges partial support from Conacyt grant No. 3255TE 
and DGAPA-UNAM grants No. IN112401 and No. IN122002. 
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